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ABSTRACT 



The "radiation inner edge" of an accretion disk is defined as the inner boundary of the region from which most of the luminosity 
emerges. Similarly, the "reflection edge" is the smallest radius capable of producing a significant X-ray reflection of the fluorescent 
iron line. For black hole accretion disks with very sub-Eddington luminosities these and all other "inner edges" locate at ISCO. Thus, 
in this case, one may rightly consider ISCO as the unique inner edge of the black hole accretion disk. However, even for moderate 
luminosities, there is no such unique inner edge as differently defined edges locate at different places. Several of them are significantly 
closer to the black hole than ISCO. The differences grow with the increasing luminosity. For nearly Eddington luminosities, they are 
so huge that the notion of the inner edge losses all practical significance. 

Key words, black holes - accretion disks - inner edge 



1. Introduction 

Accretion flows on to black holes must change character before 
matter crosses the event horizon. Two reasons account for this 
fundamental property of such flows. First, matter must cross the 
black-hole surface at th e speed of light as measured by a local in- 
ertial observer (see e.g. lGourgoulhon & Jara millo 2006), so that 
if the flow is sub-sonic far away from the black-hole (in prac- 
tice it is always the case) it will have to cross the sound barrier 
(well) before reaching the horizon. This is the property of all re- 
alistic flows independent of their angular momentum. The sonic 
surface in question can be considered as the inner edge of the 
accretion flow. 

The second reason is related to angular momentum. Far from 
the hole many (most probably most) rotating accretion flows 
adapt the Keplerian angular momentum profile. Because of the 
existence of the Inner-Most Stable Circular Orbit (ISCO) such 
flow must stop to be Keplerian there. At high accretion rates 
when pressure gradients become important the flow may ex- 
tend below the ISCO but the presence of the Inner-Most Bound 
Circular Orbit (IBCO) defines another limit for a circular flow 
(the absolute limit being given by the Circular Photon Orbit; the 
CPO). These critical circular orbits provide another possible def- 
inition of the inner edge of the flow, in this case of an accretion 
disk. 



The question is: what is the relation between the accretion 
flow edges? In the case of geometrically thin disks the sonic and 
Keplerian edges coincide and one c a n defi ne the ISCO as the 
inner edge of such disks. Paczyfiski (2000) showed rigorously 
that, independent of viscosity mechanism, presence of magnetic 
fields etc. the ISCO is the universal inner disk's edge for not 
too-high viscosities. The case of thin disks is therefore settlecQ 
However, this is not the case of non-thin accretion disks, 
i.e. the case of medium and high luminosities. The problem of 
defining the inner edge of an accretion disk is not just a formal 
exercise. Afsh ordi & Paczvnskil (12003) explored several reasons 
which made discussing the location of inner edge r = r m of 
the black hole accretion disks an interesting and important is- 
sue. One of them was, 

Theory of accretion disks is several decades old. With time 
ever more sophisticated and more diverse models of ac- 
cretion onto black holes have been introduced. However, 
when it comes to modeling disk spectra, conventional steady 
state, geometrically thin-disk models are still used, adopting 



1 In a recent paper IPenna et al.l {2010) studied the effects of mag- 
netic fields on thin accretion disk (the disk thickness Hjr < 0.07, which 
corresponds to L < 0.2 L^a). They found t hat to within a few percen t 
the magnetized disks are consistent with the lNovikov & ThorneH 19731) 
model, in which the inner edge coincides with the ISCO. 
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the classical "no torque" inner boundary condition at the 

marginally stable orbit. 
The best illustration of this fact is the case of the state-of-art 
works on measuring the black hole spin a in the microquasar 
GRS 1915+105 by fitting its observed "thermal state" spec- 
tra to these calculated (e.g. IShafee et al.ll2008t iMiddleton et alJ 
2009). These works use general relativistic version of the clas- 
sical Shakura-Sunyaev th in accretion disk model worked out by 
iNovikov & Thornel(ll973h . The Novikov-Thorne model assumes 
that the inner edge of the the disk r m = risco is also the inner- 
most boundary of the radiating region. 

Because the black hole mass of GRS 1915+105 is known 
and therefore fixed (Mo = 14 M Q + 4 M©), the surface area A of 
the radiating region, calculated in the model, depends only on 
the black hole unknown spin, a* (a* - Jc/GM 2 with / being the 
total angular momentum of the black hole). In the thermal state, 
the disk spectrum is close to that of a sum of black body con- 
tributions from different radial locations. Its shape is determined 
by the radial distribution of temperature, which in the Novikov- 
Thorne model depends on the spin, T = T(r,a*). The total ra- 
diation power L is determined by the "averaged" temperature 
To = To{a*) and the surface area A = A(a*) of the radiating re- 
gion, L = ctTqA. By calculating the spectral shape and power 
for different a* in the Novikov-Thorne model, one may find the 
best-fit estimates for the spin-dependent temperature and area. 
This is just the main idea of the spin estimate; deta i ls of the fit- 
ting are far more compl ex (see IShafee et"al]|2008t iDavis et al] 
2005; IStraub et al.ll2010|) and include, for example, a heuristic 
way of treating a contribution of scattering in accretion disks at- 
mosphere (i.e. the "hardening factor"). 

Results obtained this way bv IShafee etafl J2008h for GRS 
1915+105 showed that a* = 0.99 for the whole luminosity range 
L < 0.2 Lpm . However, for L > 0.2 L_dd, the spin estimated by 
Shafee et al. (2008) was much lower, a* » 0.8. The inconsis- 
tent spin estimates at different luminosities indicate that some 
assumptions adopted by the Novikov-Thorne model are wrong 
at high luminosities. 

This is not a surprise, because there are several physical ef- 
fects known to be important at high luminosities, but ignored 
in the classical Shakura-Sunyaev and Novikov-Thorne thin ac- 
cretion disk models. These effects are prope rly included in the 
slirn^ accretion disks models, introduced by lAbramowicz et al] 
( 1988). Advection is perhaps the best known of these "slim disk 
effects", but in the present context equally important is a sig- 
nificant stress due to the radial pressure gradient (for thin disks 
dP/dr * 0). The stress firmly holds matter well inside ISCO and 
as a result of this, at high luminosities the edge of the plunge-in 
region may be considerably closer to the black hole than ISCC0. 

Slim disks are assumed to be stationary and axially symmet- 
ric. They are described by vertically integrated Navier-Stokes 
hydrodynamical equations; no magnetic fields are considered. 
The effecti ve viscosity, believed to be generated by the MHD 
turbulence dBalbus & Hawlevlll99lh is described by the "aP" 
Shakura-Sunyaev ansatz. Figure Q] shows the slim disk location 
with respect to other analytic and semi-analytic disk models, in 



2 The names thin and slim refer to the dimensionless vertical geomet- 
rical thickness, h = H/r. For thin disks it must be h <k 1, while for slim 
disks a weaker condition h < 1 holds. 

3 Matter may be hold well inside ISCO also by magnetic 
stresses, a s pointed out by ma ny authors; see e.g. a semi-analytic 
model by [ Narayan et al. ( 2003), or MHD numerical simulations by 
iNoble. Krolik & Haw lev (2010), and references quoted in these papers. 




10-" 10- 3 0.01 0.1 1 10 10 z 10 3 10" 10 5 
optical depth at 20 M 

Fig. 1. The Figure illustrates a few best-known analytic 
and semi-analytic solutions of the stationary black hole ac- 
cretion disks. Their location in the parameter space ap- 
proximately correspond to viscosity a = 0.1 and ra- 
dius r = 20 M. For detailed reviews of these solu- 
tions see, http://www.scholarpedia.org/article/Accretion_discs 
or lKato et a l. (2008). 



the parameter space [r,h,m] described by the vertical optical 
depth r, dimensionless vertical thickness h = H/r and dimen- 
sionless accretion rate m = M/M_dd, where M_dd = 16LEdd/c 2 
is the critical accretion rate approximately corresponding to the 
Eddington luminosity (L_dd = 10 38 M/M© erg/s) in case of a disk 
around a non-rotating black holq^. 

In this paper, we discuss properties of the inner edge of slim 
accretion disks around rotating bla ck holes, using m odels sim- 
ilar to those calculated recently by ISadowskH (|2009£|. For con- 
venience, we shortly remind the slim disk basic equations in the 
Appendix A. In the following Section|2] we list six possible def- 
initions of the inner edge. These definitions reflect different (but 
partially overlapping) physical meanings and different practical 
astrophysical applications. In the following six Sections [3][8] we 
calculate the slim disk locations of these six inner edges, and dis- 
cuss their astrophysical relevance. Some of the results presented 
here have been anticipated previously by us and other authors 
in a different context of the Polish do ughnuts ( i.e. thick accre- 
tion disks; see e. g. a short review by Paczvhski 1998); see also 
iPaczvnskillfcOOOl) and lXfshordi & Paczvriskil J2003il 



2. Definitions of the inner edge 

iKrolik & Hawlevl d2002l) proposed several "empirical" defini- 
tions of the inner edge, each serving a different practical purpose 



4 Two warnings about notation, (i) Many authors use a different def- 
inition, Mem = ^Edd/c 2 . (ii) We often use the c = 1 = G convention in 
which M = r G = GM/c 2 . 

5 At http://users.camk.edu.pl/as/slimdisks a very detailed data 
base for these solutions is given. It covers the whole parameter space 
relevant for microquasars and AGN. 
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(see also the follow-up bv lBeckwith et alJ2002l) . We add to these 
a few more definitions. The list of the inner edges considered in 
this paper consists ofl 

[1] The potential spout edge r m = r pot , where the effective poten- 
tial forms a self-crossing Roche lobe, and accretion is governed 
by the Roche lobe overflow. 

[2] The sonic edge r m = r son , where the transition from subsonic 
to transonic accretion occurs. Hydrodynamical disturbances do 
not propagate upstream a supersonic flow, and therefore the sub- 
sonic part of the flow is "causally" disconnected from the super- 
sonic part. 

[3] The variability edge r m = r var , the smallest radius where or- 
bital motion of coherent spots may produce quasi periodic vari- 
ability. 

[4] The stress edge r; n = r str , the outermost radius where the 
Reynolds stress is small, and plunging matter has no dynamical 
contact with the outer accretion flow; 

[5] The radiation edge r m = r ra d, the innermost place from which 
significant luminosity emerges. 

[6] The reflection edge r m = r re f, the smallest radius capable of 
producing significant fluorescent iron line. 

In the next six Sections we discuss the six edges one by one. 
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3. The potential spout edge 

The idea of the "relativistic Roche lobe overflow" govern- 
ing accretion c lose to the black hole was first explained by 
Paczyhski (see iKozlowski et aLlll978l) . Lat er, it was explored 
in det ail by many authors analytically (e.g. lAbramowic 3 TT98TL 
1985b and by large-scale hydro dynamical simulations (e.g. 
Igumenshchev & Beloborodov 1997). It bee ame a standard con- 
cept in the black hole accretion theory. Figure [2] schematically 
illustrates the Roche lobe overflow mechanism. The left-most 
panel presents demonstrative profile of disk angular momentum 
which reaches the Keplerian value at the radius corresponding 
to the self-crossing of the equipotential surfaces presented in the 
middle panel. To flow through this "cusp" matter must have po- 
tential energy higher than the value of the potential at this point 
- such "potential barrier" is crossed only when the matter over- 
flows its Roche lobe. Precise profiles of the potential barriers and 
the angular momentum, calculated with the slim disk model, are 
presented in Figs.[3]and|4], respectively. 

The potential difference between the horizon and the spout 
is infinite, and therefore no stress may prevent the matter located 
there from plunging into the black hole. For radii greater that 
rpot, the potential barrier at r — r pot holds the matter in. Note, 
that because the dynamical equilibrium is given (approximately) 
by Vj1/ e ff = V,P/p, with p being the density, one may also say 
that it is the pressure gradient (the pressure stress) that holds the 
matter inside r pot . 

The specific angular momentum in the Novikov-Thorne 
model is assumed to be Keplerian. Slim disk models do not 
a priori assume an angular momentum distribution, but self- 
consistently calculate it from the relevant equations of hydro- 
dynamics ( IA. 1 l >- dA~8T l. These calculations reveal that the type of 
angular momentum distribution depends on whether accretion 
rate and viscosity locate the flow in the disk-like, or the Bondi- 
like type. 



6 Krolik & Hawley defined [4], [5], [6] above and in addition [7], 
the turbulence edge, where flux-freezing becomes more important than 
turbulence in determining the magnetic field structure. Magnetic fields 
are not considered for slim accretion disks, and we will not discuss [7]. 



Fig. 3. Profiles of the effective potential near the potential barrier 
(solid lines) for different accretion rates, a = 0.01 and a* = 0. 
The arrows indicate locations of the inner edge potential spout 
- r pot corresponding to the top of the potential barrier. The hori- 
zontal dashed lines present the energy of the gas overflowing the 
barrier calculated at r pot . 
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Fig. 4. Angular momentum profiles for slim disk solutions with 
a = 0.01 (left panel) and a = 0.1 (right panel). In both panels 
three curves are presented for sub-Eddingtonian, Eddingtonian 
and super-Eddingtonian accretion rates. The thin dotted line 
presents the Keplerian angular momentum profile. 



In the Bondi-type accretion flows the angular momentum is 
everywhere sub-Keplerian, X < -Lk- These flows are typical for 
high viscosities and high accretion rates, as in the case of a — 0. 1 
and m = 10 shown in Figure [4] This is the only Bondi-like 
flow in this Figure. In the disk-like accretion flows, the angu- 
lar momentum of the matter in the disk is sub-Keplerian every- 
where, except the strong-gravity region r pot < r < r cen where 
the flow is super-Keplerian, £, > The radius r cen > r^co 
corresponds to the ring of the maximal pressure in the accretion 
disk. This is also the minimum of the effective potential. The ra- 
dius r pot < nsco marks a saddle point for pressure and effective 
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r pot r isco r pot r po t 



Fig. 2. An illustrative visualisation of the Roche lobe overflow. The leftmost panel schematically presents disk angular momentum 
profile and its relation to the Keplerian distribution. The middle panel shows the equipotential surfaces. The dotted region denotes 
the volume filled with accreting fluid. The rightmost panel presents the potential barrier at the equatorial plane (z = 0) and the 
potential of the fluid (Ws) overflowing the barrier. The figure is taken from http://www. scholarpedia.org/article/Accretion_discs 



potential; this is also the location of the "potential spout inner 
edge", r m = r pot . 

Note that i n the classic solutions for spherically accretion 
flows found bv Bondil(ll952l) the viscosity is unimportant and the 
sonic point is saddle, while in the "Bondi-like" flows discussed 
here, angular momentum transport by viscosity is essentially 
important and the sonic point is usually nodal. Therefore, one 
should keep in mind that the difference between these types of 
accretion flows is also due to the relative importance of pressure 
and viscosity. For this reason a different terminology is often 
used. Instead of "disk-like" one uses the term "pressure-driven" 
and instead of "Bondi-like" one uses "v iscosity-driven" (see e.g. 
iMatsumoto etaT]|l984_[Kato et al.ll2008l) . 

From the above discussion it is clear that the location of this 
particular inner edge r pot is formally given as the smaller of the 
two roots, r ± = (r+, r_), of the equation 

[XM-XxW]^ =0. (1) 

The larger root corresponds to r cen t- Obviously, equation ([]]) has 
always a solution for the disk-like flows, and never for the Bondi- 
like flows. Figure[5]shows a division of the parameter space into 
regions occupied by Bondi-like and disk-like flows. 

The location of the potential spout inner edge r pot is shown 
in Figure [6] for a — 0.01. Note that for small accretion rates, 
m < 0.3, location of the potential spout inner edge coincides 
with ISCO. At m « 0.3, the location of the potential spout jumps 
to a new position, which is close to the radius of the innermost 
bound circular orbit, rmco- This behavior is now well-known. 
It was noticed first by Koziow skTetaLl (fl978l) for Polish dough- 
nuts, and by [Abramowi cz et alJ (Il988h for slim disks. We con- 
clude the Section on the potential spout inner edge by giving an 
approximate formula for its location, 

r pot (a*,m) = Min [ (0.275 - 0.410a* + 0.143a* 2 )»T L4 + 
4.45 - 4.87a* + 8.06a* 2 - 6.38a* 3 ; 
0.985 nsco]. (2) 

The formula (0 is valid for a = 0.01. 



4. The sonic edge 

By a series of algebraic manipulations one reduces the slim 
disk equations ( lA.lb -( TA~8l to a set of two ordinary differential 
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Fig. 5. Location of the Bondi-like and the disk-like slim accre- 
tion disks in the [a, m] parameter space. The Bondi-like accre- 
tion flows are everywhere sub-Keplerian. Disk-like flows are 
sub-Keplerian for most radii, but have also a super-Keplerian 
part around ISCO. 

equations for two dependent variables, e.g. the Mach number 
77 = — V/c 2 , = —VY.jP and the angular momentum X = -11$, 

dlnq _ M(r,i?,X) dlnL _ N 2 (r,Tt,£) 

dlnr ~~ D{r,T],£) dlnr ~ D(r,r),£) 

For a non-singular physical solution the nominators N\ and A_ 
must vanish at the same radius as the denominator T>. The de- 
nominator vanishes at the sonic edge (or sonic radius) where the 
Mach number is close to unity, i.e. 

D(r,ri,£) lr=rsm =0. (4) 

For low mass accretion rates, smaller than about 0.3MEdd in case 
of a* - 0, the sonic edge r son locates close to ISCO, indepen- 
dently on the viscosity a, as Figure|7]shows. At about 0.3MEdd a 
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Fig. 6. Location of the potential spout inner edge r pot for viscos- 
ity a = 0.01 and a* = 0. Solid lines show the exact location of 
r pot given by equation The approximation (|2]i is shown by 
dashed lines, and the location of ISCO by dotted lines. 




Fig. 7. Location of the sonic point as a function of the accretion 
rate for different values of a, for a non-rotating black hole, a* = 
0. The solid curves are for saddle type solutions while the dotted 
curves present nodal type regimes. 



plained it in terms of the disk-Bondi dichotomy. The depen- 
dence of the sonic point location on the accretion rate in the 
near-Eddington regime is more complicated and is related to the 
fact that in this range of accretion rates the transition from the 
radiatively efficient disk to the slim disk occurs near the sonic 
radius. 

The topology of the sonic point is important, because phys- 
ically acceptable solutions must be of the saddle or nodal type; 
the spiral type is forbidden. The topology may be classified by 
the eigenvalues A\ , Ai, A3 of the Jacobi matrix, 



J = 



,11) 


i)V 


(W 


dr 
(97V, 


dr, 
dN, 


f)£ 
<W, 


dr 
8N 2 


dr, 
<JN 2 


d£ 
ON, 


dr 


en 


d£ 



(5) 



Because det(JT) = 0, only two eigenvalues A\,A2 are non-zero, 
and the quadratic characteristic equation that determines them 
takes the form, 



2 A 2 - 2 A ti(J) - [tr(J" 2 ) - tr 2 (J")] = 0. 



(6) 



The nodal type is given by A1A2 > and the saddle type by 
A1A2 < 0, as marked in Figure [7] with the dotted and the solid 
lines, respectively. For the lowest values of a only the saddle 
type solutions exist. For moderate values of a (0.1 < a < 0.4) 
the topological type of the sonic point changes at least once with 
increasing accretion rate. For the highest a solutions have only 
nodal type critical points. 

The extra regularity conditions at the sonic point 
Nj(r, j], £) — are satisfied only for one particular value 
of the angular momentum at the horizon which is the eigenvalue 
of the problem. •£;« is not known a priori, and should be found. 
Figure [8] shows how does X,„ depend on the accretion rate and 
the a viscosity parameter. 




1 10 

accretion rate [M Edd ] 



qualitative change occurs, resembling a "phase transition" from 
the Shakura-Sunyaev behavior, to a very different slim-disk be- 
havior. 

For higher accretion rates the location of the sonic point sig- 
nificantly departs from ISCO. For low values of a the sonic point 
moves closer to the horizon down to ~ AM for a = 0.001. 
For a > 0.2 the sonic point moves outward with increasing 
accretion rate reaching values as high as 8M for a = 0.5 and 
lOOMndd- This effect w a s firs t noticed for small accretion rates 
by iMuchotrze b-Czernvl d 19861) and later inves t igated in a wide 
range of accretion rates bv lAbramowi cz et al.l (119881) . who ex- 



Fig. 8. Angular momentum at the horizon dependence on accre- 
tion rate for solutions with different values of a for a* = 0. 



5. The variability edge 

Axially symmetric and stationary states of slim accretion disks 
represent, obviously, only a certain theoretical idealization. Real 
disks are non-axial and non-steady. In particular, one expects 
transient coherent features at accretion disk surfaces — clumps, 
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Fig. 9. The fluid flow trajectories in slim accretion disks shown by thin solid lines for different accretion rates. Locations of r pot 
and the location of black hole horizon are shown by thick gray solid and broken lines, respectively. For small accretion rates, the 
pattern of trajectories consists of very tight spirals (almost circles) for r > r pol * r^co and very wide spirals (almost a radial fall) 
for r < rpo t . In this case, there is a sharp transition from almost circular motion to almost radial free-fall that clearly defines the 
variability edge as r var = risco- For higher accretion rates, the fluid trajectories are wide open spirals in the whole inner region of 
the flow and the variability edge makes no sense. 



flares, and vortices. Orbital motion of these features could quasi- 
periodically modulate the observed flux of radiation, mostly 
through the Doppler effect and the relativistic beaming. Let II be 
the "averaged" variability period, and All a change of the period 
during one period due to radial motion of a spot. The variability 
quality factor Q may be estimated by, 



1 _ All _ AQ _ 1 dSl _ 2^- 1 d ° M " 

q ~ ~TT ~ ~o~ ~ nd7 r ~ ^o^ "d7 i7 



(7) 



where u" '/ u l = dr/dt and u r and u' are contravariant components 
of the four velocity. The period relates to the orbital angular ve- 
locity by n = 2n/Q.. Using the relations (see Appendix A for the 
explanation of the notation used), 



V 



Va 



u 



It = 



Vl - V 2 r 

y Va Va 



rVA rVA V(l - V 2 )(l - (V^) 2 ) 



(8) 



with V being the radial velocity as measured by an observer 
corotating with the fluid, one obtains: 



1 


dlogQ 


-l 




2n 


dlogr 







(9) 



where, 

yd' 



= + [l - X - X 2 a* 2 (a 2 + 1) - X 5 a* 4 ] 



y<l> 



Vl -(V* -ojR) 2 



-1/2 



(10) 



with X = 2r c /r. From dA.2l ) and (IA.51 > it is clear than AA > 
outside the black hole horizon. Note that in Newtonian limit it 
is X <§c 1 and one has f*(a*, r) = 1. In this limit V, V* are the 
radial and azimuthal component of velocity, and the formula (0 
takes its obvious Newtonian form. 
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Fig. 10. The quality factor Q profiles for different accretion rates. 
Triangles show r pot for each rate. The vertical dashed line de- 
notes the location of ISCO. 



Behavior of the quality factor Q is shown in Figure [10] 
Profiles for four accretion rates are drawn. As Fig. [9] shows the 
lower accretion rate the smaller radial velocity component and 
therefore the quality factor Q in general increases with decreas- 
ing accretion rate. For the lowest values of m a rapid drop is 
visible at ISCO corresponding to the change in the nature of the 
flow (gas enters the free-fall region below ISCO). For higher 
accretion rates such behaviour is suppressed as the trajectories 
become wide open spirals well outside ISCO. 

Note that our definition (0 of the quality factor Q, essen- 
tially agrees with a practical definition of the variability qual- 
ity factor Qq defined by observers with the help of the observa- 
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tionally constructed Fourier variability power spectra, I(v). Here 
/(v) is the observed variability power (i.e. the square of the ob- 
served amplitude) at a particular observed variability frequency 
v. Any observed quasi periodic variability with the frequency 
~ Vo shows in the power spectrum as a local peak in I(v), cen- 
tered at a certain frequency vo. The half-width Av of the peak 
defines the variability quality factor by Qo = v/Avq. 

Quasi periodic variability with kHz frequencies, called kHz 
QPO, is observed from several low-mass neutron star and 
black hole b inarie s. In a pioneering and important research, 
iBarret et"a"D (120051) carefully measured the quality factor for 
a particular source in this class (4U 1608-52) and found that 
Qo ~ 200, i.e. that the kHz are very coherent. They argued that 
Qo ~ 200 cannot be due to kinematic effects in orbital motion of 
hot spots, clumps or other similar features located at the accre- 
tion disk surface, because these features are too quic kly sheared 
out by the differenti al rotation of the disk (see also Bath et al. 
Il974t lPringlel |l98lt). They also argued that although coherent 
vortices may survive mu ch longer times at the disk surface (e.g. 
Abramow icz et al.ll 19951) . if they participate in the inward radial 
motion, the observed variability Qo cannot be high. Our results 
shown in Figure [TOl illustrate and str engthen this point. We also 
agree with the conclusion reached bv lBarret et alJd2005l) that the 
observational evidence against orbiting clumps as a possible ex- 
planation of the phenomenon of kHz QPO, seems to point out 
that this phenomenon is most probably due to the accretion disk 
global oscillations^. F or exc ellen t reviews of the QPO oscillatory 
models see lWagoneri(ll999l) and lKatd(l200ll) . 

Although clumps, hot-spots, vortices or magnetic flares can- 
not explain the coherent kHz QPOs with Qo ~ 200, they cer- 
tainly are important in explaining the cont inuous, featureless 
Fouri er variability power spectra (see e.g. lAbrampwicz et al.l 
1199 It ISchnittmanll2005t fPechacek. Karas & Czernvl 120081 and 
references quoted there). Our results shown in Figure [TO] indi- 
cate that: (i) at low accretion rates, a sharp high-frequency cut- 
off in I(v) may be expected at about the ISCO frequency, (ii) 
at high accretion rates there should be no cut-off in I(v) at any 
frequency, (iii) the logarithmic slope p = (dlnl/dlnv) should 
depend on m. 

A more quantitat ive descrip tion of (i)-(iii) will be given in a 
future publication (IStraubll20 1 Ot) . 



6. The stress edge 

The Shakura-Sunyaev model assumes that there is no torque at 
the inner edge of the disk, which in this model coincides with 
ISCO. Slim disk model assumes that there is no torque at the 
horizon of the black hole. It makes no assumption on the torque 
at the disk inner edge, but calculations prove that the torque is 
small there. 

The zero-torque at the horizon is consistent with the small 
torque at the inner edge of slim disks, as Figure [TTI shows. The 
Figure presents the relative importance of the torque T by com- 
paring it with the "advective" flux of angular momentum Mj 
(c.f. equation IB. 11 1. For the viscosity parameter a smaller than 
about 0.01, the ratio x — T jMj both at r pot and r son is smaller 
than 0.01 even for highly super-Eddington accretion rates, and 
for small accretion rates the ratio is vanishingly small, y ~ 10~ 5 . 
For high viscosity, a - 0.5, the ratio is very small for small ac- 




1 10 
accretion rate [M Edd ] 



Fig. 11. Ratio of the angular momentum flux caused by torque 
to the flux caused by advection calculated at r pot (top) and r son 
(bottom panel) versus mass accretion rate for a number of values 
of a and a* = 0. The r pot profiles for high viscosities terminate 
when disk enters the Bondi-like regime. 



cretion rates, x < 10~ 3 and still smaller than about 0.1 even for 
super-Eddington accretion rates (calculated at the sonic radius 
as the disk enters the Bondi-like regime for such high accretion 
rates). 

To define the stress inner edge r str one has to specify the 
characteristic value of the torque parameter x- Profiles of r str 
for a few values of x and a = 0.01 are shown in Fig. [12] The 
stress edge for^ — > is located at ISCO for low accretion rates. 
When accretion rate exceeds ~ 0.3MEdd it departs from ISCO 
and moves closer to BH approaching its horizon with increasing 
m. Behaviour of r str profiles for higher (> 0.1) values of x is dif- 
ferent - they move away from the BH as the angular momentum 
profiles become flatter with increasing accretion rates (compare 
Fig.©. 

In the case of disk-like accretion with a very low viscosity 
a < 0.01, it is with high accuracy, 



'pot 



(11) 



7 IBarret et"aTJ d2005l) found also how Q varies in time for each of the 
two individual oscillations in the "twin-peak QPO". This gives strong 
observational constraint s for possible oscillat ory models of the twin 
peak kHZ QPO; see also lBoutelier et alJ d2010h . 



In this case the "inner edge" inherits both the sonic edge and 
the potential spout edge properties; suggesting a small torque. It 
looks, as this is indeed the cas e. By pushing the M HD numeri- 
cal simulations to their limits, IShafee et al.l d2008l) and recently 
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r [M] 



Fig. 12. Profiles of r str defined as the radius with given value of 
the torque parameter^ for a - 0.01. BH horizon and ISCO are 
also shown with dot-dashed and dashed lines, respectively. 

iPenna et all (l2010t) calculated a thin, H/r < 0.1, disk-like ac- 
cretion flow, and showed that for it the inner edge torque was 
small. 

7. The radiation edge 

As discussed in the previous section, the torque at r pot < osco is 
small, but non-zero and therefore there is orbital energy dissipa- 
tion also at radii smaller than ISCO. Thus, some radiation from 
this region takes place and the inner edge is not expected to co- 
incide with the radiation edge, r ra d < r pot . In Fig. [13] we present 
profiles of r lat j defined as radii limiting area emitting given frac- 
tion of disk total luminosity. For low accretion rates (< 0.3MEdd) 
disk emission terminates close to ISCO as the classical models of 
accretion disks predict. Locations of the presented r ra d are deter- 
mined by the regular Novikov & Thorne flux radial profile. For 
higher accretion rates disk becomes advective and the maximum 
of the emission moves significantly inward. As a consequence 
of the increasing rate of advection (and resulting inward shift of 
r pot ) the efficiency of accretion drops down. 

We want to stress here that the location of the radiation edge 
is not determined by the location of the stress edge (as some au- 
thors seem to believe), but by the fact that significant advection 
flux brings energy into the region well below ISCO. 

Let r out » rc be the outer radius of the disk. The total lumi- 
nosity of the disk could be estimated from 

L = Me mi + (QQ.) mA -Me oM -(Q£l) oM 

= M£ md + Q md - M£ out - gout- (12) 

It is £2 out w 0, e out ~ 0, and from this one derives 

L = M [e rad + xUn\nd] = W (13) 

where x is the ratio of the viscous torque to the advective flux of 
angular momentum (see Fi gures fTTI and [TZt . 



Because ^ <K 1, the efficiency of accretion r\ depends mainly 
on the specific energy at the inner edge, e, al j. The further away 
is the inner edge from ISCO (and closer to the black hole), the 
smaller is the efficiency. 




2 4 6 8 10 



r [M] 

Fig. 13. "Luminosity edges" defining inner radii of area emitting 
given amount of the total disk radiation. The lines are drawn 
for 95%, 99% and 99.9% of the total emission. The dashed line 
shows the location of the potential spout inner edge r pot . The 
gravitational suppression of the radiation has been taken into ac- 
count. 



8. The reflection edge 

The iron Kq, fluorescent line is one of the characteristic fea- 
tures observed in many sources with black h ole accretion disks 
dMmedl2006tlRemillard & McClintockl2006l) . The intensity and 
the shape of this line depends strongly on the physical condi- 
tions close to t he inner edge. This was di scussed by many au- 
thors, including lRevnold s & Fabian (2008) who gave three con- 
ditions for line formation: (i) the flow has to be Thomson-thick 
in the vertical direction; (ii) disk has to be irradiated by exter- 
nal source of X-rays (hard X-ray irradiation plays crucial role in 
the process of fluorescence and changes the ionization degree of 
matter); (iii) the ionization state should be sufficiently low (iron 
cannot be fully ionized). 

We point out here, that this first condition is sufficient for the 
formation of the reflection continuum, but formation of the fluo- 
rescent iron line requires an even stronger condition i.e. that the 
effective optical depth of the flow should be higher than unity: 

Teff = V T abs(T a b.s + T es ) > 1. (14) 

This is because, the fluorescence requires efficient absorption of 
high energy photons by iron ions. In Fig.[l4]we present profiles 
of the effective optical depth T e g in different regimes of accre- 
tion rates for a = 0.1 and a* - 0. Three characteristic types of 
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their behaviour are shown: sharp drop, maximum and monotonic 
at the top, middle and bottom panels, respecively. Behaviour 
for different values of a and a* is qualitatively similar (but not 
quantitatively as in general r e ff increases with decreasing a). 
Top panel, corresponding to the lowest accretion rates, shows 
a sharp drop in r e ff near ISCO. The same beh avior was noticed 
previously e.g. by iRevnolds & Fabian] (120081) . The drop could 
clearly define the inner reflection edge r re f ~ risco limiting the 
radii where formation of the fluorescent iron line is prominent. 
The middle panel, corresponding to moderate accretion rates, 
shows a maximum in r e ff near ISCO. The non-monotonic be- 
haviour is caused by the fact that regions of moderate radii out- 
side ISCO become radiation pressure and scattering dominated. 
Note, that the top of the maximum of r e ff stays near ISCO in 
a range of accretion rates, but for accretion rates greater than 
0.3MEdd it moves closer to the black hole with increasing m as 
the disk emission profile changes due to advection. The bottom- 
most panel corresponds to super-Eddington accretion rates. The 
profiles are monotonic in T e s and define no characteristic inner 
reflection edge. Close to the black hole such flows are effectively 
optically thin reaching T e g — 1 on about few tens of gravitational 
radii. 

When effective optical depth of the flow becomes less then 
unity, our approximation of radiative transfer by diffusion with 
grey opacities (Eq. IA.8b becomes not valid. In such case full ra- 
diative trans fer through accretion disks atmospheres sho uld be 
solved (e.g. iDavis et all 120051: iRdzahska & Madeill2008l) . Still, 
our results allow us to estimate roughly how far from the black 
hole the iron line formation is most prominent, assuming that 
disk is uniformly illuminated by an exterior X-ray source. For 
accretion rates smaller than 0.3MEdd, the reflection edge is lo- 
cated very close to ISCO and we may identify shape of the iron 
line with gravitational and dynamical effects connected to ISCO. 
In case of higher but sub-Eddington accretion rates, the maxi- 
mum of the effective optical depth is located inside ISCO what 
may possibly allow us to study extreme gravitational effects 
on the iron line profile. However, the assumption that the line 
is formed at ISCO is no longer satisfied. The super-Eddington 
flows have smooth and monotonic profiles of the effective opti- 
cal depth. Therefore, the reflection edge cannot be uniquely de- 
fined and no relation between shape of the fluorescent lines and 
ISCO exists. Finally, one should keep in mind that such lines can 
be successfully modeled by clumpy absorbing mate rial and may 
have nothing to do with relativistic effects (see e.g. Mille r et al.l 
2009, and references therein). The role of the ISCO in deter- 
mining the shape of the Fe lines was also questioned in the past 
(based on different reasoning) by IRevn olds & Begelman ( 1997) 
whose arguments were then refuted bv lYoung et al.l (Il998l) . 



9. Conclusions 

We addressed the inner edge issue by discussing behavior of six 
differently defined "inner edges" of slim accretion disks around 
the Kerr black hole. We found that the slim disk inner edges 
behave very differently than the corresponding Shakura-Sunyaev 
and Novikov-Thorne ones. The differences are qualitative. Even 
for moderate luminosities, M > 0.3 Mem, there is no unique 
inner edge. Differently defined edges locate at different places. 
For nearly Eddington luminosities, the differences are huge and 
the notion of the inner edge losses all practical significance. 

We summarize the properties and locations of the six inner 
edges in Table Q] It refers to a* = 0, but the qualitative behavior 
is similar for a* + 0. 



We conclude, by showing in Figure [TBI differences between 
the Shakura-Sunyaev and slim-disk (in the disk-like case) treat- 
ment of the inner disk physics. The innermost part of a Shakura- 
Sunyaev disk is shown in the left column in Figure Q3] and the 
innermost part of a slim disk is shown in the right column. The 
upper panel shows angular momentum in the disk (the solid 
line) in reference to the Keplerian distribution (the dashed line). 
ISCO, indicated by the dash-dotted line is at the radius where 
the Keplerian angular momentum has its minimum. The poten- 
tial spout (a square) and the center (a triangle) are defined as 
crossings of the angular momentum in the disk line with the 
Keplerian line. For slim disks they are at two different radii, on 
both sides of the ISCO. For Shakura-Sunyaev disks they merge 
into one singular location at ISCO. The lower panel shows the 
cross section of the disk. The slim disk has everywhere a finite 
thickness while the Shakura-Sunyaen disk is singular at ISCO (it 
has a zero thickness there). The sonic radius (a cross) is where 
the accretion component of the velocity equals the local sound 
speed. In slim disks, the sonic point corresponds to a critical 
point of the set of differential equations, that through the regu- 
larity conditions defines the global eigensolution of the problem. 
The Shakura-Sunyaev disk is described by local algebraic equa- 
tions and this global eigenvalue aspect is missing, thus location 
of a sonic point is of no relevance. 
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Appendix A: The Kerr geometry slim disks 

The Shakura-Sunyaev models are local: they are described by 
algebraic equations, valid at any particular (radial) location in 
the disk, independently of physical conditions at different loca- 
tions. Contrary to that, the slim disk models of accretion disks 
are non-local. They are described by differential equations glob- 
ally connecting physical conditions at all radial locations from 
infinity to the black hole horizon. 



Initially, mode ls of slim disks have been constructed by 
lAbramowicz et all (1 19881). who used the pseudo-Newtonian po- 
tential of Paczvriski & Wiita (1 19801) and Newto nian equations 
derived by Paczvriski & Bisnovat vi-Koganl(ll98ll) and later im- 
proved by Muchotrzeb & Paczvriski (1982), Matsumoto et al. 
(1 19841) and iMuchotrzebl dl983l) . General relativistic version 
(the Kerr metric) of the slim d i sk equations was derive d 
and elaborated by iLasotal (1 19941) . lAbramowicz et al.1 d 19961) 
Gammie & Popham (1 19981) . and most recently by ISad owski 
(2009) who made several corrections and improvements to re- 
sults of the previous authors, and who numerically constructed 
slim disk solutions in a wide range of parameters applicable 
to the X-ray binaries. In particular, he calculated the solutions 
in the whole relevant range of accretion rates, from very sub- 
Eddingtonian, to moderately super-Eddingtonian on es. In this 
paper we follow notation and conventions adopted by Sadowski 
( 2009). The Kerr geometry slim disk equations adopted here are: 



(i) The mass conservation: 



M 



-2ttZA 



l/2_ 



V 



Vi -v 2 



(A.l) 



r+h 

where Z = J pdz ~ 2Hp is disk surface density and V is the 
gas radial velocity as measured by an observer at fixed r who 
co-rotates with the fluid. Here 



A = r -IMr-a 1 . 



(A.2) 



(For the Kerr metric description see e.g. iKato et all I2008L or 
any textbook on general relativity). Equation dA.U has the same 
form in the Shakura-Sunyaev model. 

(ii) The radial momentum conservation: 



V dV \_dP_ 
\-V 2 ~dr ~ ~r ~ Ydr 



where 



MA (O-Q+XCl-Q,-) 



1 - Q 2 R 2 



r 4 + r 2 a 2 + 2M ra 2 , 



(A.3) 



(A.4) 



(A.5) 



Q. = u^/u' is the angular velocity with respect to the stationary 
observer, Q = Q - a> is the angular velocity with respect to the 
inertial observer, = +M'' 2 /(r 3 ' 2 ±aM 1 ^ 2 ) are the angularfre- 
quencies of the co-rotating and counter-rotating Keplerian orbits 
and R = A/(r 2 A 1 ^ 2 ) is the radius of gyration. In the Shakura- 
Sunyaev model this equation is a trivial identity = because 
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the radial pressure and velocity gradients vanish, and rotation is 
Keplerian, Q. — Q.^ . 

(iii) The angular momentum conservation: 



M A 1/2 A 1/2 y 

—U-£ in ) = -aP 

In r 



(A.6) 



where X. = -u<p is the specific angular momentum, y is the 
Lorentz factor and P - 2Hp can be considered to be vertically 
integrated pressure. The const ant a is the standard alph a vis- 
cosity parameter introduced by Shakur a & Sunvaevl(ll973l) . The 
constant X/« is the angular momentum at the horizon, unknown 
a priori. It provides an eigenvalue linked to the unique eigenso- 
lution to the set of slim disk differential equations constrained by 
proper boundary and regularity conditions. The algebraic equa- 
tion ( IA.6b is the same as in the Shakura-Sunyaev model, except 
that the Shakura-Sunyaev model assumes that -£;« = Xi(ISCO). 
(iv) The vertical equilibrium: 



S// 2 



t l - a\e l - 1) 
2r4 



(A.7) 



with e = u, being the conserved energy associated with the time 
symmetry. The same equation is valid for the Shakura-Sunyaev 
model. 

(v) The energy conservation: 



aPAy 1 dQ. 32 oT 4 



M 



1 



dr 



dp 

3 kI. 2nrp Tj, - 1 \ dr 



T,^) (A.8) 



here T is the disk central temperature. The right hand side of 
this equation represents the advective cooling and vanishes in 
the Shakura-Sunyaev model. Because in the Shakura-Sunyaev 
model rotation is Keplerian, Q. = Q£, which means that Q. is 
a known function of r and therefore the first term on the left- 
hand side (which represents viscous heating) is algebraic. The 
second term, which represents the radiative cooling (in diffusive 
approximation) is similar in the Shakura-Sunyaev model. 

Appendix B: No torque at the black hole horizon 

The assumption about (vanishingly) small torque in the region 
between black hole and accretion disk is well motivated phys- 
ically. Let us recall that the very meaning of a torque Q is 
that it transports angular momentum without transporting mass. 
Correspondingly, the total angular momentum flux j through a 
surface equals, in general, 



J = Mj + Q, 



(B.l) 



where M is the mass flux, and j is the angular momentum density 
(per unit mass). However, the torque is only a phenomenological 
concept. Microscopically, the flux j should be seen as a differ- 
ence of material fluxes that come from the opposite sides of the 
surface, j = [M + j + ] - [M_ _/'_]. One also has M - M + - M+, and 
j = (j+ + jJ)/(j + M+ + 7_M_). Microscopically then, the torque is 
equal Q = 2M+ M-(j+ - j-)/(M+ + M_). It necessarily vanishes 
when all matter crosses the surface in only one direction, i.e. 
when either M + — or M_ = 0. As the only one-side matter flux 
is the fundamental property of the black horizon, one concludes 
that there sho uld be no torque at th e black hole surface. 

Since the Blandford & Znaiek] d 19771) process energizes the 
jet (and disk) by extracting rotational energy of a black hole 
through a kind of electromagnetic braking, some astrophysicists 
argue that in this case there must be a "Maxwell" torque be- 
tween the black hole and outside matter. However, by looking at 



the Blandford-Znajek process from the quantum electrodynam- 
ics perspective, one sees only ingoing, but not outgoing photons. 
Thus, there is only one-way traffic of photons, and no torque pos- 
sible. The photons with negative energy and angular momentum 
that are present in the ergosphere, are responsible for the slow- 
ing down the hole, similarly to negative energy particles in the 
classic Penrose process that must necessarily have also a nega- 
tive angular momentum. This point of view, that the Blandford- 
Znajek process is an electromagnetic version of the Penrose pro- 
cess, was recently discussed in context of the classical Maxwell 
electrodynamics (in Kerr geometry) by seve ral authors, in par- 
ticular most forcefully bv lKomissarovl d2008l) . 

Here, we generalize Komissarov's point to any material field, 
not only the electromagnetic one. Following Komissarov, let us 
consider the local ZAMO (or FIDO) observer in the Kerr geom- 
etry. His four velocity in terms of the Killing vectors rf (time 
symmetry) and (axial symmetry) is given by n' = q(jf + cl>%'), 
where u is the angular velocity of frame dragging, and q > 
follows from normalization n'n k gik = — 1. Let us now consider a 
general matter or field, described by an unspecified stress-energy 
tensor T' k . The energy flux in the ZAMO frame is E' = -T' k n k . 
The energy acquitted by the black hole is 



I 



T'n" dNi > 0, 



(B.2) 



where J dN; is the surface integral over the horizon. The in- 
equality sign follows from the fact that the locally measured en- 
ergy must be positive. The above integral may by transformed 
into 



< E 



j qT'ttf +cot k )dN i = q H {E a 



(B.3) 



where the index H denotes horizon, and Eco and J m are the 
"energy at infinity" and the "angular momentum at infinity" ac- 
quired by the black hole absorbing the corresponding fluxes of 
these quantities defined by, 



(B.4) 



From (IB. 31 one concludes that E^ > J^ojh- As in the classic 
Penrose process, the necessary condition for the extraction of 
energy at infinity is that the energy (at infinity) absorbed by a 
black hole is negative, £«, < 0. This is equivalent to JcoOJh < 0. 
Thus, in a way fully analogous to the line of arguments that is 
made discussing the Penrose process, one may say that if energy 
at infinity increases because the black hole absorbed a negative- 
at-infinity energy, then the black hole must also slow down by 
absorbing matter or electromagnetic flux with negative angular 
momentum. 




Fig. 14. Profiles of the effective optical depth r e ff for a — 0.1 and 
a = in three different regimes of accretion rate. Vertical lines 
denote locations of the BH horizon (dotted) and ISCO (dashed). 
Three types of behavior of T e g are seen: sharp drop at ISCO 
for smallest accretion rates, maximum near ISCO for moderate 
accretion rates, and monotonic everywhere for highest accretion 
rates. 
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